Test of many-electron QED effects in the hyperfine splitting of heavy high-Z ions 
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A rigorous evaluation of the two-photon exchange corrections to the hyperfine structure in lithiumlike heavy 
ions is presented. As a result, the theoretical accuracy of the specific difference between the hyperfine splitting 
values of H- and Li-like Bi ions is significantly improved. This opens a possibility for the stringent test of the 
many-electron QED effects on a few percent level in the strongest electromagnetic field presently available in 
experiments. 
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Accurate measurements of the ground-state hyperfine structure performed in H-like 209 Bi, 165 Ho, 185 Re, 187 Re, 207 Pb, 203 T1, 
and 205 T1 GH1 were intended to probe QED in the strong electromagnetic field generated by a heavy nucleus. However, 
accurate calculations revealed that the uncertainty of the predicted hyperfine splittings, which mainly originates from the nuclear 
magnetization distribution correction (the Bohr-Weisskopf effect), is comparable in magnitude with the QED correction, see 
e.g. Refs. |0]. Accordingly, a direct identification of QED effects on the hyperfine splitting in heavy H-like ions appeared to 
be unfeasible. It was shown instead, that this uncertainty can be significantly reduced in a specific difference of the hyperfine 
splitting values of H- and Li-like ions with the same nucleus 0]: A'E = AE (2s ^ -£AE( ls \ where AE^ and AE^ are the 
hyperfine splittings of H- and Li-like ions, respectively, and the parameter £ is chosen to cancel the Bohr-Weisskopf correction. 
The parameter £ can be calculated to a rather high accuracy independently of the employed nuclear magnetization distribution 
model. Thereby, the stringent tests of QED in strong fields can be achieved by studying the specific difference of the hyperfine 
splitting values in H- and Li-like ions. 

Till recently there existed only an indirect measurement of the 2s hyperfine splitting in lithiumlike Bi ion with an accuracy of 
about 3% [8]. Direct measurements with high-precision laser spectroscopy are feasible at the current experimental storage ring 
(ESR) and future HITRAP facilities in GSI [9]. Just recently, after 13 years of various attempts, the hyperfine splitting of the 
ground state Li-like Bi ion has been directly observed in GSI lioll . 

Achievement of the required theoretical accuracy for the specific hyperfine splitting difference for H- and Li-like heavy ions 
demands the rigorous evaluation of various QED and interelectronic-interaction effects. Since the influence of one-electron 
QED corrections is considerably reduced in the specific difference, the total value of A'E is essentially determined by the 
screened radiative and interelectronic-interaction corrections. Recently, the screened self-energy and a domina nt p art of the 
screened vacuum-polarization contributions have been calculated rigorously within a systematic QED approach 111 ill 1211 . This 
calculation represents an essential advance beyond the local screening potential approximation employed in the previous works 
lfl3l - [l5ll . As concerns the interelectronic-interaction contribution, up to now it was evaluated rigorously only up to the first order 
in 1/Z, see, e.g., Ref. [16]. The contributions of second- and higher-order in 1 /Z were calculated within the Breit approximation 
employing many-body perturbation theory and configuration-interaction methods I13lll5lll7lll8l1 . however, for high-Z ions such 
calculations can provide only an approximate estimation of these corrections. In the present Letter we report on the complete 
evaluation of the second-order interelectronic-interaction corrections within a rigorous QED approach. As the most interesting 
application of these results we present improved theoretical predictions for the specific difference between the ground-state 
hyperfine splitting values of H- and Li-like Bi ions. 

The second-order interelectronic-interaction corrections in the presence of an external potential correspond to the third-order 
perturbation theory terms. Nowadays, several approaches are used for derivation of the formal expressions for perturbation se- 
rial terms from the first principles of QED: the two-time Green-function method Ill9ll . the covariant-evolution-operator method 
02oll . and the line profile approach lEHl . Here, we employ the two-time Green-function method. To simplify the derivation 
we specify the formalism regarding the closed shell electrons as belonging to a redefined vacuum. In this way we have to 
consider all two-loop diagrams for the valence electron in the presence of magnetic perturbation, i.e. 30 nonequivalent dia- 
grams. These diagrams encompass the second-order interelectronic-interaction corrections, the one-electron two-loop, and the 
screened one-loop radiative corrections. In our recent works 111 ill 1211 we have evaluated the screened radiative corrections. The 
generic types of the second-order interelectronic-interaction diagrams, where we now focus on, are depicted in Figs. [T] [2] The 
interelectronic-interaction corrections to the hyperfine splitting were usually denoted by the terms B(aZ) JZ and C(Z, aZ) /Z 2 , 
see e.g. Ref. 01511 . The term B(aZ)/Z determines the interelectronic-interaction correction of the first order in 1/Z, while the 
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FIG. 1 : Feynman diagrams representing the three-electron part of the two-photon-exchange corrections to the hyperfine splitting. The wavy 
line indicates the photon propagator and the double line indicates the electron propagators in the Coulomb field. The dashed line terminated 
with the triangle denotes the hyperfine interaction. 
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FIG. 2: Feynman diagrams representing the two-electron part of the two-photon-exchange corrections to the hyperfine splitting. Notations are 
the same as in Fig.Q] 



interelectronic-interaction corrections of second- and higher-orders are denoted by the term C(Z, aZ) /Z 2 . Here, we isolate the 
terms of the second- and higher-orders: C(Z, aZ)/Z 2 = C(aZ)/Z 2 + D(Z, aZ)/Z :i , where C(aZ)/Z 2 corresponds to the 
interelectronic-interaction corrections of the second-order in 1/Z only (Figs. [UI9, an d the term D(Z,aZ)/Z 3 represents the 
third- and higher-order corrections in 1/Z. As in Ref. [ 19], we call the parts of diagrams reducible when an intermediate-state 
energy coincides with the reference-state energy and irreducible otherwise. The irreducible parts of the three-electron diagrams 
depicted in Fig. Q~]yield the following contributions 
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where the prime on the sums over intermediate states indicates that terms with vanishing denominators should be omitted in the 
summation. The irreducible parts of the two-electron diagrams depicted in Fig.|2]yield 
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where the prime on the sums indicates that in the summation we omit the reducible and infrared-divergent terms, namely, those 
with e ni + e ri2 = e a + e b in the ladder-W diagrams, with e ni — Ep a , £n 2 — £ Qb in the direct parts of the cross-W diagrams, 
with £ ni = e„ 2 — e a ,e b in the exchange parts of the cross-W diagrams, with e ni + £„ 2 = E a + £ b and s ni + £„ 3 = £ a + £ b and 
£« 2 = £ n 3 = £Qb - Ap a Q a in the ladder-S diagrams, with e ni = e Pa , £„ 2 = e Qb and e ni = e Pa , £„ 3 = e Qb and £„ 2 = £„ 3 = 
£gb in the direct parts of the cross-S diagrams, with e ni — £„ 2 = e a ,s b and £ m = £„ 3 = e a , £ b and £„ 2 = £„ 3 = e a ,£ b in the 
exchange parts of the cross-S diagrams. In Eqs. (HJ-®, a and b refer to the valence- and core-electron states, respectively; the 
sum over b runs over all closed-shell states, P and Q are the permutation operators giving rise to the signs (— l) p and (— 1)*^ of 
the permutation, respectively. I(uS) is the interelectronic-interaction operator [ 19], u = 1 — iO preserves the proper treatment of 
poles of the electron propagators. The energy difference A ni „ 2 is defined as A„ lJl2 = e ni — e„ 2 . To is the electronic part of 
the hyperfine-interaction operator and G a is the multiplicative factor depending on the quantum numbers of the valence electron 
(see for details Ref. [ 15]). The modified wavefunction |£) is defined as follows 

IO = £° . (9) 

£a £n 

n 

We refer to all contributions, where the energies of intermediate states and the reference state coincide, the infrared-divergent 
contributions, and the nondiagrammatic terms as the remaining reducible part. Formal expressions for this part, which are rather 
bulky, will be published elsewhere. Following the treatment of Ref. B22I1 we introduce a nonzero photon mass. In such a way we 
regularize the infrared divergences and cancel them analytically. 

Now let us discuss the numerical evaluation procedure. The infinite threefold summations over the spectrum of the Dirac 
equation have been performed employing the dual-kinetic -balance finite basis set method 12311 with basis functions constructed 
from B-splines [24]. The Fermi model for the nuclear charge density and the sphere model for the magnetic moment distribution 
have been utilized. The summations over magnetic substates have been performed analytically by means of standard formulas 
and also numerically as an independent check. The most problematic part has consisted in evaluation of the two-electron terms, 
which contain the integration over the energy of the virtual photon oj. In order to avoid strong oscillations arising for large 
real values of uj, we have performed a Wick rotation of the integration contour, as it was done in the calculations of two- 



photon exchange corrections to the Lamb shift Refs. I 25H27I1 . and have employed the integration contours such as in Ref. B27I1 . 
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TABLE I: Individual contributions to the two-photon exchange correction C(aZ)/Z 2 for the ground-state hyperfine structure of the Li-like 

209gj80+ 



Contr. 


Feynman 


Coulomb 


3el, A 


0.001685 


0.002229 


3el, B 


-0.001942 


-0.002489 


3el, C 


0.001154 


0.001036 


3el,D 


0.003781 


0.003914 


2el, lad-W 


0.003401 


0.003960 


2el, cr-W 


0.000363 


-0.000019 


2el, lad-S 


0.001155 


0.001226 


2el, cr-S 


0.000207 


-0.000001 


reducible 


-0.009063 


-0.009116 


Total 


0.000740 


0.000740 



TABLE II: Individual contributions to the specific difference A' E for 209 Bi in meV. 





AE (2s) 


£A£ ,(ls) 


A'E 


Dirac value 


844.829 


876.638 


-31.809 


Interelectronic interaction 








~ l/Z 
~ l/Z 2 


-29.995 




-29.995 


0.258 




0.258 


~ l/Z 3 and higher orders 


-0.003(3) 




-0.003(3) 


QED 


-5.052 


-5.088 


0.036 


Screened QED 


0.193(2) 




0.193(2) 


Total 






-61.320(4)(5) 



However, special care should be taken for the identification of the pole structure of the integrands, because they are essentially 
more complicated than in the case of the two-photon exchange corrections to the energy levels. 

In what follows we present our result for the case of Li-like Bi utilizing the following values for the nuclear properties: 
(r 2 ) 1 / 2 = 5.5211 fm Q, F = 9/2-, and fx = 4.1106(2)//jv H- We have performed calculations in both Feynman and 
Coulomb gauges and the corresponding individual contributions to the C(aZ)/Z 2 are presented in Table J] The total result is 
gauge independent on the level of the numerical accuracy, that serves as an accurate check for both the derived formulas and 
the numerical procedures. As an additional check, we have reproduced the third-order many-body perturbation theory (MBPT) 
result, considering the Breit approximation in the derived expressions. Finally, we have found that the result of the rigorous QED 
evaluation of the two-photon exchange correction 0.000740 is in reasonable agreement with the l/Z 2 term 0.00075 extracted 
from the large-scale configuration-interaction Dirac-Fock-Sturm calculation. 

Now let us come to the consideration of the specific difference between the ground state hyperfine splitting values for H- 
like and Li-like Bi. The cancellation of the Bohr-Weisskopf effect appears with £ = 0.16886 for the case of Bi. In Table El 
we present the current status of individual contributions to the specific difference. The rigorous evaluation of the two-photon 
exchange corrections improves the accuracy of the interelectronic -interaction term by an order of magnitude in comparison with 
previous calculations 111 ll The accuracy of the screened QED contribution has been also increased due to recent rigorous 
evaluations of the Wichmann-Kroll parts of the electric and magnetic loops [30], which was accounted for in Refs. Ill U ll 2 1 
within some approximations. Thus, the remaining theoretical uncertainty for the specific difference comes from the uncalculated 
Wichmann-Kroll parts of the screened vacuum-polarization correction and from the l/Z 3 and higher orders interelectronic - 
interaction term. The second uncertainty in the total value of A'E arises from the uncertainty of the nuclear magnetic moment, 
the nuclear polarization corrections Il32ll . and other nuclear effects, which do not completely cancel in the specific difference. 
It should be noted that the nuclear magnetic moment uncertainty can be larger due to a chemical shift 13 ill . Employing the 
experimental value of the Is hyperfine splitting A-Eexp = 5.0840(8) eV id and the theoretical result for the specific difference, 
one can easily find the hyperfine splitting for Li-like Bi Ai?( 2s ' = 797.16(14) meV, where the accuracy is fully determined 
by the uncertainty of the experimental value. As one can see from Table [Til the one-electron QED corrections are strongly 
canceled in the specific difference and the dominant QED contributions comes from the many-electron effects. Therefore, the 
theoretical accuracy achieved now for the specific difference allows to test the many-electron QED effects on the few percent 
level, provided the hyperfine splittings in H- and Li-like bismuth are measured with a relative accuracy of about 10~ 6 . When 
the QED corrections will be tested and found to be valid, the comparison between the theoretical and experimental values will 
enable the determination of the nuclear magnetic moments and their volume distribution. 

In summary, we have rigorously calculated the two-photon exchange correction to the hyperfine splitting in heavy Li-like ions. 
As a result we have significantly increased the accuracy of the specific difference, thus providing the theoretical prerequisite for 
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a test of many-electron QED effects at strongest electromagnetic fields. Further extensions of these calculations to the g factor 
of Li-like and B-like heavy ions may also serve for an independent determination of the fine structure constant from QED at 
strong fields II 3 311 . 
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